Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 1 February 2008 



(MN MfeK style file vl.4) 



The X— ray Spectrum of the Dwarf Nova SS Cyg in 
quiescence and outburst 

C. Done 1 and J. P. Osborne 2 

1 Department of Physics, University of Durham, South Road, Durham, DH1 3LE 

2 Department of Physics & Astronomy, University of Leicester, University Road, Leicester. LEI 7RH. 
C*** ' 1 February 2008 

o\ ■ 



ABSTRACT 



We reanalyse archival GINGA and ASCA X-ray data from SS Cyg in both qui- 
escence and outburst, using multi-temperature plasma models for the continuum and 
line emission, together with their reflection from the X-ray illuminated white dwarf 
and accretion disk. Reflection is clearly detected in all the spectra, but its contribution 
is larger in the softer X-ray spectra seen in outburst than in quiescence. This sup- 

■ ports models in which the quiescent inner disk is not present or not optically thick, 
' so that the only reflector is the white dwarf surface rather than the white dwarf plus 

fT) , the disk. The amount of reflection in outburst is also more consistent with the hard 

■ X-rays forming a corona over the white dwarf surface rather than just an equatorial 
C*~) ' band. We detect partially ionized absorption in the ASCA outburst spectrum, which 

is probably the X-ray signature of the outflowing wind. The ASCA data also allow a 
detailed measure of the elemental abundances. We find that all detectable lines from 
the hot plasma (except perhaps Si and S) are a factor ~ 2.5 x weaker than expected 
from solar abundances. We examine possible deviations from coronal equilibrium, but 
conclude that the heavy elements are truly underabundant in SS Cyg. 

In quiescence the underlying intrinsic spectrum is consistent with a single temper- 
^ ■ ature plasma. This conflicts with the expected cooling of the material, and requires 

^ | that either the plasma is reheated or that the cooler emission is masked by absorp- 

■ tion in an optically thin inner disk which affects only that part of the boundary layer 

emission below the disk. Such absorption is also required if the data are to match the 
theoretical models (Narayan & Popham 1993) for the quiescent hard X-ray emission. 

. By contrast, the outburst spectrum is much softer and is dominated by the cooling 

' components, so that it cannot be fit by a single temperature model. There are no read- 

ily available (or believable) theoretical models for the hard X-ray emission in outburst, 
but the emission observed can be described by material with a continuous temperature 
distribution. Although the outburst spectrum is much softer than the quiescent spec- 
trum, the total luminosity of the X-ray component in both cases is similar, showing 
that the optically thin plasma emission is important even in outburst. 

Key words: stars: individual: SS Cyg - cataclysmic variables - binaries: close - 
X-rays: accretion - X-rays: accretion 



1 INTRODUCTION 

The non-magnetic Cataclysmic Variable stars are systems 
where a white dwarf accretes via a disk from a Roche lobe 
filling, low mass, companion star. Dwarf novae (DN) are a 
subclass of these systems which show optical/UV outbursts 
of 2-6 magnitudes which are thought to be triggered by a 
disk instability giving a transient increase in the accretion 
rate (see e.g. the reviews by Cordova 1995, Warner 1995 
and Osaki 1996). As the accreting material falls inwards, 



the transition from the the disk velocity field to that of the 
white dwarf will lead to up to half of the accretion energy 
being dissipated (Pringle 1977, Kluzniak 1987). 

Detailed theoretical models of the boundary layer emis- 
sion are not well developed due to the inherent complexity of 
modelling the strongly shearing, probably turbulent transi- 
tion region between the disk and white dwarf surface. How- 
ever, general characteristics of the emission can be derived. 
Material leaves the inner edge of the disk with a tempera- 
ture ~ 10 5 K, and then is heated (presumably by shocks) 
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as it joins the boundary layer. If the boundary layer is op- 
tically thick then the gas cools efficiently and so collapses 
quickly onto the white dwarf surface, forming a geometri- 
cally thin equatorial belt around the white dwarf. The lumi- 
nosity emerges at the local blackbody temperature, which 
is of order 20 — 50 eV for the case where all the remain- 
ing accretion luminosity is radiated (Pringle 1977; see also 
Frank King and Raine 1992). Conversely, if the gas is opti- 
cally thin then radiative cooling is much less efficient, the 
gas can be heated up to much higher temperatures before 
reaching equilibrium with radiative cooling and/or adiabatic 
expansion, giving temperatures of order ~ 10 keV. Since the 
optical depth of the accreting gas is strongly linked to the 
mass accretion rate this predicts that mainly EUVE/soft X— 
rays are expected from DN in outburst, while hard X-rays 
can be produced by DN in quiescence (Pringle and Savonijc 
1979; Tylenda 1981; King and Shaviv 1985, Narayan and 
Popham 1993; Popham and Narayan 1995). A logical exten- 
sion of this model is that the vertical density gradients in the 
disk could ensure that some fraction of the accreting mate- 
rial remains optically thin even in outburst, producing some 
hard X-ray emission, as observed (Patterson and Raymond 
1985a; 1985b). 

The resulting hard X-ray spectrum in both quiescence 
and outburst should certainly consist of multi-temperature 
components as the accreting gas must cool in order to set- 
tle onto the white dwarf photosphere. Thus the continuum 
should be characterised by continuous temperature distribu- 
tion, with the contribution from each temperature compo- 
nent being weighted by its emissivity. However, in addition 
to this intrinsic hard X-ray continuum, there should also be 
a component formed from reflection of this spectrum from 
the illuminated white dwarf surface and/or accretion disk. 
The reflection albedo is energy dependent, as at low energies 
the scattering probability is decreased by photo-electric ab- 
sorption, whereas at higher energies Compton downscatter- 
ing and the reduction of the scattering cross-section deplete 
the number of photons reflected. This results in a broad 
spectral "bump" , centered around 20 keV, while the depen- 
dence on photo-electric absorption at low energies gives a 
pronounced feature at the iron K edge energy in the re- 
flected spectrum and an associated Fe Ka fluorescence line, 
both of which are a function of the ionization state and el- 
emental abundances of the reflecting gas (e.g. Lightman & 
White 1988; George & Fabian 1991; Matt, Perola & Piro 
1991). Such models have been applied to magnetic CV sys- 
tems (Beardmore et al 1995; Done, Osborne and Beardmore 
1995; Cropper, Ramsey and Wu 1996; van Teeseling, Kaas- 
tra and Heise 1996), but have not yet been extended to the 
non-magnetic systems. 

A further possible distortion of the hard X-ray spec- 
trum in DN is from complex absorption. P Cygni line pro- 
files are seen in these systems in the UV SilV, CIV and 
NV resonance lines, but only in outburst, when the accre- 
tion rate is high (see e.g. the recent review by Drew 1993). 
The mass loss rate and ionization structure in this wind is 
poorly known, but the models which have most success in 
matching the blueshifted UV absorption features are opti- 
cally thick at soft X-ray energies (e.g. Mauche and Raymond 
1987; Hoare and Drew 1991; 1993). Thus partially ionized 
absorption may play an important role in determining the 
observed soft X-ray spectrum in outbursting DN spectra, 



though the wind declines rapidly after the outburst peak 
(Verbunt et al 1984), so this is unlikely to be important for 
the quiescent emission. 

These spectral distortions will bias the results from 
simple thermal model fits to the data. This is important 
as recent exciting developments in accretion disk structure 
have attempted to include a self consistent description of 
the boundary layer emission. In the quiescent state the disk 
density is low, so local cooling may be inefficient enough 
that radial energy transport is also important. Numerical 
solutions by Narayan and Popham (1993) including these 
advected energy terms give predictions of the temperature- 
cmissivity structure of the boundary layer. However, as both 
reflection and ionized absorption give an observed spectrum 
that is flatter than the intrinsic continuum, their neglect 
can give an overestimate of the amount of high temperature 
material. 

SS Cyg is the brightest DN at in hard X-rays, so is 
the obvious candidate for detailed spectral studies. Here we 
analyse the available archival data from GINGA (2 spec- 
tra) and ASCA (1 spectrum). These show the system in 
three different states, quiescence, decline from a normal (fast 
rise) outburst, and in anomalous (slow rise) outburst, re- 
spectively. We evaluate the importance of reflection and 
partially ionized absorption, and so obtain the underlying 
intrinsic continuum temperature-emissivity distribution of 
the hot plasma. These data have been previously reported 
by Yoshida et al (1992), Ponman et al (1995) and Nousek 
et al (1994) but these authors only fit the spectra with sim- 
ple thermal models. Here we include reflection and ionized 
absorption, and explicitally fit continuous temperature dis- 
tribution models to the underlying intrinsic spectral shape. 

We clearly detect the signatures of reflection in all the 
spectra. In outburst this is consistent with the hard X-ray 
emitting plasma forming a corona over the white dwarf, 
forming a reflection spectrum from both the white dwarf 
surface and inner disk. However, in quiescence there is sig- 
nificantly less reflection, which we interpret as indicating 
that the inner disk is disrupted or optically thin during qui- 
escence. The ASCA data also show partially ionized absorp- 
tion, as expected from the outflowing wind. 

The abundances of the elements can be constrained, es- 
pecially in the ASCA data. The line emission is significantly 
weaker than expected from solar abundances. We examine 
possible deviations from coronal equilibrium, but conclude 
that none of these can cause the observed line deficit. The 
high iron abundance relative to the other elements in the 
Anders & Grevasse (1989) 'solar' compilation is inconsis- 
tent with the data. Newer 'solar' abundance determinations 
(Grevasse, Noels & Sauval 1996) give a relative iron abun- 
dance much closer to the Morrison & McCammon (1983) 
values, which are similar to the elemental abundance ratios 
in the ASCA data presented here. 

The underlying continuum is well described by a con- 
tinuous (power law) temperature-emissivity distribution of 
hot material. In quiescence this is strongly weighted towards 
high temperatures, so that it is indistinguishable from a sin- 
gle temperature component. This is inconsistent with the 
expected cooling of hot gas. Either there is some reheating 
process or the cool components are masked by partial ab- 
sorption, perhaps by an inner disk which is optically thin to 
electron scattering which photoelectrically absorbs the soft 
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X-ray emission from the boundary layer below the disk. By 
contrast, the outburst spectrum is dominated by cool com- 
ponents, and cannot be described by a single temperature 
plasma. While the emission is significantly softer, the total 
luminosity is similar to that in quiescence, assuming that 
the material cools down onto the photosphere. 



2 THEORETICAL MODELS 
2.1 Continuum Emission 

If gas is heated up to some temperature T max and then 
allowed to cool down to the photospheric temperature T bb 
at constant pressure in a constant gravitational potential so 
that the luminosity is solely driven by the temperature drop 
then the resulting spectrum can be described by a sum of 
thermal spectra from T ma x to Tbb such that 

where e(T, n 2 ,v) is the bremsstrahlung spectral emissivity 
and e(T, n 2 ) is the total (integrated over v) bremsstrahlung 
emissivity. More general cooling flow models include a 
weighting factor factor (T/T max ) a in the integrand (see e.g. 
Johnson et al 1992; Done et al 1995). 

Explicit calculations of the emission from the boundary 
layer are extremely difficult due to the strongly shearing and 
probably turbulent nature of the flow. Any accretion disk 
problem is inherently two dimensional (at least) as the mat- 
ter flow is radial while the energy is transported vertically 
(through radiation) as well as radially (through advection). 
Reducing the problem to 1 dimension by vertical averag- 
ing of the accretion disk structure (the standard Shakura- 
Sunyaev approach) makes the equations far more tractable. 
This is the approach used by Narayan & Popham (1993) 
and Popham & Narayan (1995) to calculate the temperature 
structure of the boundary layer. Their models then give the 
predicted temperature-emissivity distribution in the bound- 
ary layer. Figures la and b show spectra calculated from 
the models of Narayan & Popham (1993) for the quiescent 
emission for a white dwarf of mass M = Mq and radius 
R = 5 x 10 8 cm for a mass accretion rate of M = 10 -9 ' 5 
and 10 -10 ' 5 Mq yr" 1 , respectively. The dotted lines in these 
figures show a comparison with an a = cooling flow model 
of maximum temperature equal to that of the highest tem- 
perature emission present in the advective boundary layer 
model, which is 8.78 and 16.4 keV, respectively. The two 
models are remarkably similar at low M, but begin to di- 
verge at higher mass accretion rates. In outburst however, 
their models give no hot optically thin component to the 
emission (Popham & Narayan 1995), probably because of 
the limitations of the 1 dimensional approach whereas the 
hard X-ray emission is thought to be from material escap- 
ing from the boundary layer to form a corona over the white 
dwarf/disk. 

Calculations in (at least) 2 dimensions are required in 
order to follow coronal structure. Kley (1989) and Hujeirat 
(1995) solve numerically a fully 2 dimensional set of (simpli- 
fied) equations describing the hydrodynamic and thermal 
properties of the gas. The results of Kley (1989) are es- 
pecially relevant as these include most of the white dwarf 
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Figure 1. Spectra calculated from the models of Narayan & 
Popham (1993) for the quiescent emission from a IMq object 
accreting at logM = —9.5 (a: top panel) and — 10.5Mq yr -1 (b: 
lower panel) are shown by the solid lines. The dotted lines give the 
comparison spectrum expected from plasma of solar abundance 
heated to the same maximum temperature as in the Narayan & 
Popham (1993) models, cooling by radiation alone, at constant 
pressure in a constant gravitational potential. The models are 
very similar, especially at low mass accretion rates. 

surface in the calculation, and so can follow any hot, opti- 
cally thin gas which expands to form a corona around the 
white dwarf, while the grid used by Hujeirat (1995) only 
covers the equatorial belt of the white dwarf. Kley (1989) 
does indeed show that an optically thin, hot (~ 6 keV for 
his assumed system parameters) corona forms around the 
white dwarf and inner disk during outburst. The tempera- 
ture and density contours can in principle be combined to 
give the temperature-emissivity distribution (which is the 
observed quantity) but in practice the published data plots 
are too complex to convert into this form so we cannot di- 
rectly compare these calculations with the data. 

However, even a two dimensional approach is probably 
not adequate, as this cannot describe oblique shocks in the 
plane of the disk which are perhaps the most likely out- 
come of the disk/star interaction. Even more seriously, all 
these models (both 1 and 2 dimensional) assume that the 
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viscosity in the boundary layer and the disk is produced 
the same mechanism (generally assumed to be turbulent 
convective motions) and so use the same parameterisation 
to describe the viscosity in both regions. This is at odds 
with recent results which show that convective turbulence 
transports angular momentum inwards not outwards as is 
required for accretion (Balbus, Hawley & Stone 1996), and 
that the magnetic dynamo models (Tout & Pringle 1992) 
based on the magnetohydrodynamic instabilities (Balbus & 
Hawley 1991) are the most likely source of the strong viscos- 
ity inferred for accretion disks. The viscosity in the boundary 
layer may then be very different from that in the disk since 
it is not at all obvious that the MHD dynamo can work out- 
side of the ordered velocity field in the disk (since it uses 
Keplarian shear velocities to transform radial magnetic field 
to azimuthal, the Parker instability to transform azimuthal 
field into vertical, and the Balbus-Hawley instability to re- 
generate radial field from vertical. The dissipation required 
for viscosity then comes from magnetic field reconnection 
from opposing vertical field lines: Armitage, Livio & Pringle 
1996). 

2.2 Reflection from the White Dwarf Surface 

The intrinsic continuum can be reflected by any optically 
thick material. The reflection models calculated to date have 
generally assumed the specific geometry of a (pointlike or 
extended) source illuminating a flat slab (Lightman and 
White 1988; Matt, Perola and Piro 1991; George and Fabian 
1991; Zycki et al 1994; Magdiarziz and Zdziarski 1995; Van 
Teesling et al 1996). Here the geometry of the boundary 
layer emission is more complex, as we expect that the X- 
ray emission will be produced in an equatorial belt around 
the white dwarf. The vertical extent of this belt is model de- 
pendent. The advective boundary layer models of Narayan 
and Popham (1993) have h/r ~ 0.3 — 1 (R. Popham, private 
communication), while King and Shaviv (1984) and Pringle 
and Savonije (1979) envisage the quiescent boundary layer 
as extending as a corona over the whole of the white dwarf 
surface. In outburst all the models agree that the optically 
thick boundary layer should be extremely compact both in 
height and in radial extent. However, the structure of the 
residual optically thin hard X-ray emission that can be seen 
in outburst is less clear, and such high temperature gas is 
again likely to expand adiabatically out of the equatorial 
plane before it cools (Pringle and Savonije 1979; King and 
Shaviv 1984). Thus we use the vertical extent of the bound- 
ary layer as a free parameter. 

The radial extent of the boundary layer is observed to 
be very small, < 1.15Rwd, from eclipse mapping of the X- 
ray emission in the quiescent dwarf novae HT Cas (Mukai 
et al 1996). Interestingly, the width of the radiative bound- 
ary layer in the models of Narayan and Popham (1993) is 
1.15Rwd for mass accretion rates of M = 1CP 9,5 Mq yr _1 , 
and increases for smaller mass accretion rates. This is in po- 
tential conflict with the data as the hard X-ray flux from 
HT Cas only requires an accretion rate of ~ 2 x 10~ 12 M Q 
yr _1 (Mukai et al 1996), but the fractional contribution of 
the hard X-rays to the bolometric luminosity is highly un- 
certain and could be very small (Van Teeseling, Beuermann 
and Verbunt 1996). Here we follow the observations and as- 
sume that the boundary layer is thin in the radial direction. 



The boundary layer emission will clearly illuminate the 
white dwarf surface, but the extent to which this can be 
seen depends on the accretion disk geometry. In outburst 
it is generally supposed that the disk is optically thick and 
has an inner radius which is very close to the white dwarf 
surface (e.g. Popham and Narayan 1995). However, there 
is considerable debate on the structure of the accretion disk 
during quiescence. Several models have suggested that there 
is no inner disk in quiescence, in order to explain the delay in 
the onset of the UV flux rise during outburst with respect 
to the optical (as seen by e.g. Mauche et al 1995). Meyer 
and Meyer-Hofmeister (1994) proposed that the inner disk 
could be evaporated by a coronal siphon flow; Livio and 
Pringle (1992) proposed that the inner disk formation during 
periods of low mass accretion rates was inhibited by low 
level (B ~ 10 4 — 10 5 G) magnetic fields, while Armitage 
et al. (1996) proposed that the inner disk is depleted by a 
magnetic wind. Alternatively, the inner disk may be present 
in quiescence, but not be optically thick. This is suggested 
by Mukai et al (1996) to explain the marginally narrower 
X-ray than optical eclipse ingress/egress width in HT Cas. 
If the accretion disk has a column of ~ 10 23 - 24 C m~ 2 then it 
is optically thick to photoelectric absorption of the 2-3 keV 
flux (which dominates the ASCA countrate) emitted by the 
boundary layer below the disk but is optically thin to the 
optical emission from the white dwarf photosphere as long 
as the disk density is less than ~ 10 15-16 cm~ 3 so that the 
free-free opacity is small. Thus we assume that an optically 
thick disk extends down to the boundary layer in outburst, 
but that in quiescence the inner disk may be disrupted or 
optically thin. 

The proposed geometries are shown in Figure 2a and b. 
The boundary layer half thickness subtends an angle (3 at the 
center of the white dwarf, and the system is viewed at some 
inclination angle i (measured from the normal to the disk 
plane). The appendix shows that the reflection component 
from the white dwarf surface in SS Cyg, where i = 37° 
(Cowley et al 1980), can be approximated by that from an 
isotropically illuminated slab viewed at 60°, irrespective of 

2.3 Reflection from the Disk 

There should be an additional reflection spectrum from the 
optically thick disk if it does extend down to the inner 
boundary layer (see figure 2a). This is not so simple to quan- 
tify, firstly as the geometry is complex (the reflection spec- 
trum from the region of the disk shadowed by the white 
dwarf cannot be seen by the observer) and secondly because 
the amount of reflection depends also on the distribution of 
incident photon angles, and is larger for large incident an- 
gles such as are predominantly expected here. Photons which 
skim along the surface of the disk have a large path length 
through the material. They are are then absorbed/reflected 
in the upper layers of the disk, so it is very easy for any 
reflected photons to escape (e.g. George and Fabian 1991). 
There is also the problem that the reflector intercepts not 
just the intrinsic spectrum but also the reflection spectrum 
from the white dwarf surface. Thus we will just treat the lim- 
iting cases of an infinitely thin disk illuminated by a bound- 
ary layer with /3 = and (3 = 7r/2. The former case is similar 
to that of a point source above a disk, as for each point on 
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Figure 2. The assumed geometry of the boundary layer. The X— 
ray emission forms a corona over the white dwarf surface, extend- 
ing up to an angle /3 towards the pole. The observer is assumed 
to be on the Y axis, viewing the system at some inclination i with 
respect to the normal to the accretion disk in the X-Y plane. In 'a' 
this disk is assumed to be optically thick all the way down to the 
boundary layer, whereas in 'b' the inner disk is either not present 
or is optically thin so that the emission from the boundary layer 
below the equatorial plane can also be seen 

the boundary layer near the equator the emission is limited 
to a hemisphere pointing away from the white dwarf sur- 
face, and half of this will intercept the disk while the other 
half will not. Thus the effective solid angle subtended by 
the disk is Q/2-k = 1. In the latter case, where the emission 
extends all the way up to the poles, the reflection spectrum 
normalisation is smaller. Radiation emitted at the poles is 
again confined to a hemisphere pointing away from the white 
dwarf surface but this does not intersect the disk at all, 
and so contributes only to the continuum emission rather 
than to the disk reflection. In fact the solid angle subtended 
by the disk to an optically thick X-ray emitting sphere is 
S1/27T = 1/2 (Chen and Halpern 1989), so as (3 -> tt/2 the 
additional reflection component is then 0.5 x that expected 
from an isotropically illuminated disk viewed at inclination 
i. For SS Cyg the angles are such that an 0/27T = 1/2 disk 
reflection component viewed at 37° is almost identical (in 
shape and normalisation) to that from the white dwarf sur- 
face {Q./2-k = 1) viewed at the mean angle of 60°, so the 
additional reflection spectrum expected from the inner disk 
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(if it exists) is 1 — 2 x that expected from a slab viewed at 
60°. 



3 OBSERVATIONS 

SS Cyg was observed twice with GINGA (2-20 keV), once 
in quiescence (25-27 Nov 1987) and once during decline 
of a normal outburst (2-5 Dec 1990). The data were ex- 
tracted in 64 second binning using standard selection crite- 
ria (maximum pointing offset from the source of 0.4°, soft 
electron rates of SOL < 10, magnetic rigidity between 7 < 
COR< 20, surplus count rate above the upper discrimina- 
tor of SUD< 10, and Earth elevation angle between 6 < 
YELEV< 120). 

The GINGA proportional counters are split into two 
separate layers, top and mid, where the mid layer is more 
sensitive above 15 keV, but has very little effective area be- 
low 8 keV (Turner et al 1989). The data mode chosen for 
the majority of both SS Cyg observations was MPC2, in 
which the top and mid layer are combined, rather than the 
usual faint source data mode of MPC1 where the mid layer 
spectra can be discarded to increase signal-to-noise. The 
standard 'universal' background subtraction technique, in 
which blank sky observations over a period of several months 
around the date of the observation are used to calculate the 
contribution from the diffuse cosmic X-ray background, par- 
ticle background and the induced satellite radioactive decay 
components (Hayashida et al 1989; Williams et al 1992) is 
optimised for MPC1 mode data, so is not necessarily appro- 
priate here. Also, SS Cyg is in the galactic plane (b = —7) 
so the spectrum at the lowest energies can also be contam- 
inated by the diffuse galactic background. Instead we use a 
'local' background (Williams et al 1992), where data from a 
nearby (in both time and galactic coordinates) were also ex- 
tracted, using the same selection criteria as the source except 
that a larger range in geomagnetic parameters (SUD< 14) 
and offset angle (< 0.5°) is allowed. This blank sky data 
is then used to model the background components in the 
source data. For the 1987 data the background observation 
is not sufficiently long to sample all the same background 
levels as seen during the source observation so we exclude 
parts of the SS Cyg observation with SUD> 8. Two nearby 
backgrounds were used to model the 1990 data as the back- 
ground taken between the two main blocks of source data 
was noisy. Systematic errors induced by using only a short 
segment of background data are included by increasing the 
error bar on the background by a factor 1.5. This is a con- 
servative procedure and probably results in the overly small 
values of \ 2 derived in Sections 4.1 and 4.2. The background 
subtracted data is then attitude corrected in order to pro- 
duce the final dataset, and a systematic error of 0.5 per cent 
is added to take account of residual uncertainties in the de- 
termination of the instrument response. 

The default (Rev 1) screened ASCA (0.6-10 keV) data 
for each instrument were taken from the archive. These were 
extracted with elevation angle > 10°, geomagnetic rigidity 
> 6 GeV/c, pointing within 0.01° of the target position, not 
in the South Atlantic Anomaly or in fine guidance mode. 
For the SIS data there is the further criteria that the bright 
Earth angle is > 20° , and all faint mode data are converted 
to bright mode, and combined with the original bright mode 
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data. Source spectra for the SIS and GIS were extracted over 
a circle of radius 4 and 6 arcminutes, respectively, and back- 
ground was taken from the same image as the source data. 
The effective area curve for each instrument was derived 
from the source position data and this, together with the the 
standard response matrices (s0clg0234p40el_512v0_8i.rmf, 
Slc3g0234p40el_512v0_8i.rmf, gis2v4_0.rmf, gis3v4_0.rmf) 
give the instrument response. The source data were then 
grouped to a minimum of 20 counts per bin so that \ 2 fit- 
ting is appropriate. 



4 SPECTRAL FITTING 

The data were fitted using the XSPEC spectral fitting pack- 
age version 9.00 (Arnaud 1996). All error bars or upper lim- 
its quoted are A\ 2 = 2.7 (90 % confidence level for 1 free 
parameter). We caution that some parameters are strongly 
coupled (especially in the low spectral resolution GINGA 
data), so the allowed strengths of e.g. the Gaussian line 
and reflection continuum are strongly (anti)correlated. The 
equivalent widths of any additional Gaussian line compo- 
nents are quoted by re-evaluating the model on a dummy 
response matrix with 1000 energy grid points between 0.2 
and 20 keV. This means that the 6.4 keV line is calculated 
with respect to the line-free continuum at 6.4 keV, whereas 
with the instrumental energy binning the 'continuum' can be 
contaminated by the hot plasma line emission. The interstel- 
lar column density to SS Cyg is fixed at Nh — 3 x 10 19 cm~ 2 , 
as determined from the depth of narrow UV absorption lines 
seen against the bright outburst spectrum (Mauche, Ray- 
mond and Cordova 1988). 

For single temperature hot plasma spectra we use the 
mekal model in XSPEC (see the XSPEC users manual for 
a complete description of the supported models), while for 
a power law emissivity-temperature distribution plasma we 
use a modification of the cemekl code in XSPEC, where the 
luminosity of each temperature component is weighted by 
a factor (T/T max ) a • The code as released calculates the in- 
tegral over temperature logarithmically, without modifying 
the step size from dT to dlogT. Thus the a produced by 
the code is 1 larger than the true weighting factor. In all 
that follows the a quoted is the XSPEC weighting factor. 
Numerically the code also had to be modified by changing 
the order of integration (starting with the highest tempera- 
ture component rather than the lowest) so that the spectrum 
changes smoothly with changing T max . The density in both 
plasma codes is fixed at n — 1 cm~ 3 . Reflection of the con- 
tinuum spectrum is calculated from the code of Magdziarz 
and Zdziarski (1995), with the inclination angle fixed at 60° 
as appropriate for SS Cyg (see section 2), but with the frac- 
tional contribution left as a free parameter. All these models 
assume that the abundances with respect to H of all the el- 
ements scale together, and that this fractional abundance 
relative to the Solar determination of Anders and Gervasse 
(1989) is a free parameter. The iron fluorescence line that 
should accompany the reflection continuum is included as a 
separate narrow Gaussian, with line energy of 6.4 keV be- 
ing appropriate for a mainly neutral reflector. These four 
model components (single temperature plasma, power law 
temperature-emissivity distribution plasma, reflection con- 



tinuum and Gaussian emission line) are denoted by T, plT, 
R and g respectively in all the tables. 

There is also the possibility of complex absorption. In 
quiescence the inner disk may be present but be optically 
thin to electron scattering (section 2.2). This may give rise 
to partial absorption of the X-ray emission. The fraction of 
the emission that is occulted by the inner disk can be de- 
rived from the difference in the emitting area observed with 
and without a disk (i.e. the difference in normalisation fac- 
tors in equations 1 and 2 of the appendix). For SS Cyg at 
an inclination of 37° this gives rise to a maximum expected 
covering fraction of 0.5 if the boundary layer is concentrated 
in the equatorial plane ((3 — > 0) or a minimum covering frac- 
tion of 0.2 for P — > 90°. We assume that the optically thin 
disk is neutral and model this by the pcf code in XSPEC. 
The partially ionized absorption expected in outburst from 
transmission of the X-ray spectrum through the outflowing 
wind is modeled by the absori code in XSPEC. The ion 
populations are found by balancing the temperature depen- 
dent dielectronic and radiative recombination rates against 
photoionization as in Done et al (1992). The thermal bal- 
ance equations are not calculated self consistently. Instead 
the temperature of the material is fixed at 5 x 10 4 K, as ap- 
propriate for the derived low ionization states (see section 
4.3 and e.g. Kallman and McCray 1982 figures 21 and 22). 
The code assumes that the ionization balance is determined 
by a power law X-ray spectrum. The continuum from SS 
Cyg can be roughly characterised by a power law of pho- 
ton index P = 2 in outburst and Y — 1.8 in quiescence, so 
we fix the index at the relevant value for each spectrum. 
Model calculations of the ionization state of the wind have 
tended to concentrate on photoionization by the putative 
boundary layer blackbody emission but the harder X-ray 
emission will also contribute to the ionization balance. This 
especially important in determining the abundances of the 
ions which absorb in the X-ray bandpass, such as OVII. 
Both absorption models use the fixed abundances of Morri- 
son & McCammon (1983) and are denoted by pcf and pia 
respectively in all the tables. 

4.1 GINGA 2-20 keV Spectrum (1987): 
Quiescence 

Results from all the spectral fits in this section are detailed 
in table 1. Since these data are the same as those used by 
Yoshida et al (1992) we first check for consistency by us- 
ing their model of a bremsstrahlung spectrum with narrow 
Gaussian emission line. This gives a good fit to the data 
(xl = 15.3/22), with derived continuum and line parame- 
ters which are consistent with those quoted by Yoshida et 
al (1992). However, a bremsstrahlung spectrum is physically 
unlikely as there will be elements other than just pure hydro- 
gen in the accreting material. Thus we replace this ad hoc 
continuum and line with a hot plasma spectrum in which 
the line emission is calculated self consistently, and fix the 
absorption at the interstellar value. This gives a poor fit 
to the spectrum, with xt = 39.6/24 for a temperature of 
15.2 ± 0.7 keV and abundance of 0.71 ± 0.10. Including a 
second temperature component, as has often been suggested 
(Yoshida et al 1992; Nousek et al 1994) does not give a signif- 
icantly better fit (xl = 32.6/22 for fcTi = 9.2±f% keV and 
kT-2 = 57^37 keV). Adding further separate temperature 



© 0000 RAS, MNRAS 000, 000-000 



The X-ray Spectrum of SS Cyg 7 



components does not result in any significant improvement 
in the fit. 

Such ad hoc multi-temperature models can approxi- 
mate the spectrum of a continuous temperature distribu- 
tion. A power law temperature-emissivity distribution hot 
plasma model gives a fit which is not significantly differ- 
ent (xt — 37.6/23) from the single temperature description, 
with the data requiring that the most weight is given to the 
higher temperature components (a > 1.6). 

As the theoretical continuum spectral shape is not 
known a priori we are justified in trying an arbitrary emis- 
sivity/temperature relation. One such model describes the 
differential emission measure, 4>(T), as an exponential func- 
tion of an feth order Chebyshev polynomial, Pu (T) , such that 
4>(T) oc exp(T,akPk{T) (Lemen et al. 1989), where the expo- 
nential function ensures positivity. We modified the cp6mkl 
model in XSPEC so that the maximum temperature was 
a free parameter rather than being fixed at 10 keV, rede- 
fined the Chebyshev polynomial coefficients to a more stan- 
dard form (e.g. Press et al 1988), and again reversed the 
order of integration. Fitting with this model resulted in a 
xt = 33.0/18, similar to that from the two temperature 
plasma description. The derived emissivity distribution has 
two peaks at similar temperatures to those found in the two 
temperature model so we do not tabulate the results of this 
fit. 

All the above plasma models give systematic residuals 
at the iron line and edge energy. These are shown in figure 3 
for the best fitting two temperature plasma model, and are 
very suggestive of the presence of a reflection continuum and 
its associated iron Ka line emission. Including an additional 
cold iron line to even the single temperature plasma model 
gives a highly significant reduction in xt to 15.0/23. This 
line should also be accompanied by a reflection continuum, 
but the strength and shape of this continuum is dependent 
on the elemental abundances and ionization state of the re- 
flector as well as the geometry of the source with respect to 
the reflector and the inclination of the reflector with respect 
to the observer (see Section 2). We fix the inclination at 60°, 
and constrain the abundances in the reflector to be equal to 
those of the hot plasma. This gives only a marginal improve- 
ment in the fit, with xt = 13.2/22 compared to the 6.4 keV 
line fit described above, but physically it gives a coherent 
picture of the source as the amount of 6.4 keV fluorescent 
line is then consistent with a reflection origin. 

Replacing the single temperature continuum model 
with a power law temperature-emissivity distribution and 
fitting this, together with its reflection spectrum and a cold 
iron line as above, gives no significant improvement in the fit 
(xt = 12.2/21), with the incident spectrum being weighted 
towards the high temperature components (a > 1.5). The 
total intrinsic 25eV-100 keV flux in this model is 1.2 x 10~ 1() 
ergs cm~ 2 s~\ which gives a bolometric luminosity of the 
optically thin gas of 1.4 x 10 32 ergs s _1 assuming a distance 
of 100 pc (Bailey 1981). This implies a minimum mass accre- 
tion rate of M = 10~ 10 ' 4 Mq yr _1 assuming that the mass of 
the white dwarf is O.8M (Cowley et al 1980) and its radius 
is Rwd — 7 x 10 s cm (Nauenberg 1972). The temperature- 
emissivity distribution derived from the data is inconsistent 
with the a — 1 gas cooling model and the advective bound- 
ary layer model at these accretion rates (see Figure 1), while 
freezing a at unity gives significant excess emission at low 



o 



a 



2 5 10 

channel energy (keV) 



20 



to a two temperature plasma model. There is an excess at ~ 6.4 
keV, together with a decrement at 7-8 keV. These features are 
consistent with the iron fluorescent line and edge expected from 
reflection in cold material. 



energies. The mass accretion rate given above may be under- 
estimated by its neglect of possible optically thick emission 
from the boundary layer, but such corrections are likely to 
be small as the optically thick emission in quiescence seems 
to be weak (Ponman et al 1995). One possible way for the 
spectra to be made consistent with these models is to add in 
excess absorption. The data clearly are already well fit statis- 
tically but the addition of partial covering by cold material 
with the absorbed fraction fixed at 0.2 (see Section 3) gives 
an excellent fit (xt = 10.2/20) for a column of N H ~ 2 x 10 22 
cm~ 2 with the derived temperature-emissivity relation fixed 
at a = 1. 

Alternatively, the excess absorption can be described by 
complete covering by partially ionized material (the absori 
model in XSPEC with the power law index fixed at 1.8). 
Again, this gives a good fit (xt = 10.3/20) with the a — 1 
continuum model for a column of Nh ~ 5 x 10 21 cm~ 2 . 
Both absorption models increase the intrinsic (absorption 
corrected) 25eV-100 keV flux to 1.4 x 10 10 ergs cm" 2 s" 1 . 
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Table 1. Quiescent (GINGA) data 



model 



kJ a (keV) 



EW (eV) 



Xt 



IT 
2T 

plT 
IT+g 
lT+R+g 
plT+R+g 
pcf( P lT+R+g)0 
pia(plT+R+g) 



15.2 ±0.7 
9 2+ 29 

a -3.0 
57 +43 

-37 
iq 7+5.7 
iy -'-4.8 

15.6 ±0.7 

i3.ol 3 2 ;° 

15.4 4 



4+8.8 
-5.4 
-11 
-5.7 
21+ 11 



2 V 



0.71 ± 0.10 

-1+0.14 
-0.09 



0.69 1 



0.74 ±0.11 
0.42 ± 0.13 

0.35 4 

0.35^ 

0.38^ 

0.38^ 



-+0.15 
-0.10 
-+0.I6 
-0.11 

-,+0.16 
-0.11 
j+0.17 
-0.11 



4.1 x 10" 2 
2.7 x 10- 2 
1.7 x 10~ 2 

0.20 
4.3 x 10" 2 

4.2 x 10~ 2 
0.21 
0.10 
0.10 



2.8 



+00 
1.2 



2 9^°" 



170 ± 60 



5+ - 7 


125+ 95 
i ^°-90 




6 +0 - 9 


wo+11 




6+ 10 
u -°-0.6 


85+r 


16 +1 7 7 


6+ 11 







18 



O e 

+ 1000 
18 



39.6/24 
32.6/22 

37.6/23 
15.0/23 
13.2/22 
12.2/21 
10.2/21 
10.3/20 



a This is the maximum temperature in the plT models 

b The normalisation of the plasma model in units of 10 _14 /(47rD 2 ) X Emission Measure 
c Units of 10 21 cm" 2 

d Defined as L/nr 2 where r is distance from material of density n illuminated by a source of luminosity L 
e Covering fraction fixed at 0.2 
f Parameter fixed 



4.2 GINGA 2-20 keV Spectrum (1990): Decline 
from Normal Outburst 

These data were reported by Ponman et al (1995), who noted 
that the spectrum is significantly fainter and softer than that 
seen in quiescence. They split the observation into two sec- 
tions in order to compare them with simultaneous ROSAT 
survey data. However, they noted that there was no sig- 
nificant change in hard X-ray spectral shape between the 
two segments, even though the softer ROSAT data, which 
are dominated by a separate low temperature component, 
show substantial variability. Thus here we co-add all the 
data to get maximum signal to noise. A single tempera- 
ture bremsstrahlung continuum and (broad) Gaussian emis- 
sion line gave results consistent with those of Ponman et al 
(1995). 

We test the same series of models as described above, 
and details of all the spectral fits are given in table 2. Again, 
neither single nor multiple temperature models give a good 
fit to the continuum. However, unlike the quiescent data, 
adding a 6.4 keV iron emission line to the single tempera- 
ture continuum does not produce a statistically acceptable 
fit to the data. Similarly, a single temperature continuum 
together with a reflection continuum and line fails to fit 
the data (xt = 51.1/22). A good phenomenological descrip- 
tion of the data can only be achieved by adding a 6.4 keV 
emission line to an intrinsically multi-temperature contin- 
uum, such as a two temperature model (xl = 19.0/21), 
or a continuous power law temperature-emissivity distribu- 
tion (xt = 18.4/22). The lower value of the temperature 
weighting index a in this latter fit indicates the much softer 
spectrum seen in outburst compared to that in quiescence. 

Since the continuous temperature distribution seems 
more physically likely as well as giving a better fit to the 
data we use this as our intrinsic spectrum and include its 
corresponding reflection continuum in the fit as well as the 
6.4 keV line emission. This gives a significant reduction in xl 
to 12.2/21, so that in these data the reflection continuum it- 
self is detected separately from the line emission. The model 
here assumes that the reflecting material is neutral, so that 
the reflected edge is at 7.1 keV. This is at first sight incom- 



patible with the results of Ponman et al (1995), where the 
fitted edge energy of 8.3 keV implied highly ionized mate- 
rial. However, with the limited resolution of GINGA, the 
derived energy of spectral features is strongly dependent on 
the assumed continuum form (see e.g. Beardmore et al 1995 
for a similar example of this effect), 

Including a warm absorber (with illuminating power law 
index fixed at F = 2) does not result in a significant improve- 
ment to the fit (xt = 10.2/19). We do not tabulate the re- 
sults of this fit as the parameters of the partially ionized 
column are largely unconstrained. The intrinsic (absorption 
corrected) 25eV-100 keV flux is 6.8 x 10" 10 ergs cm" 2 s" 1 
and 5.2 x 10" 11 ergs cm -2 s _1 for models with and without 
ionized absorption, respectively. 



4.3 ASCA 0.6-10 keV Spectrum (1990): 
Anomalous Outburst 

We fit all 4 ASCA datasets simultaneously, and full details 
of the derived parameters are given in table 3. The overall 
shape is similar to that of the GINGA decline from out- 
burst data, i.e. much softer and fainter than the quiescent 
spectrum in the 2-10 keV band. As shown by Nousek et al 
(1994), neither a single temperature plasma nor two plasma 
components give a good description of the spectrum, with 
xl = 3799/1159 and 1331/1157, respectively. A three tem- 
perature model is required in order to give a statistically 
acceptable fit (xt = 1145/1155). A continuous (power law) 
temperature-emission measure distribution gives an even 
better fit for fewer free parameters, with xl = 1117/1158. 
This result is contrary to the assertion of Nousek et al (1994) 
that the lack of intermediate temperature plasma lines indi- 
cates that the emission measure - temperature distribution 
is not monotonic. The apparent discrepancy is probably due 
to their their assumption of solar abundances which led them 
to expect large line equivalent widths from any intermediate 
temperature components, while the line-to-continuum ra- 
tios in the data are a factor 2 weaker than predicted for solar 
abundance models. Their preferred phenomenological model 
fit to the data is a three temperature plasma, in which two 
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Table 2. Decline from normal outburst (GINGA) data 



model 



kT a (keV) 



A 



N b 



R 



EW (eV) 



IT 
2T 

plT 
IT+g 
lT+R+j 
2T+g 

plT+g 
plT+R+ 



6.9 ±0.3 0.71 ±0.08 



4.6 ±0.6 

loots 

7 1+0.3 
'■ L -QA 

6 4+ ' 9 

1 - lo -0.64 



1 0+0.18 



0.94 ±0.13 
0.42 ±0.12 
0.41 ±0.11 
0.41 H 



fO.21 
-0.14 



9.6 4 



0.501°;!* 

35+° 13 



1.34 x 10" 2 
1.13 x 10~ 2 
3.18 x lO- 3 
3.05 x 10" 2 
1.39 x 10" 2 
1.39 x 10~ 2 
5.2 x 10~ 3 
1.2 x 10~ 2 
2.21 x 10" 2 
2.34 x 10" 2 



1 00+ 54 



^1+ - 35 
U ' Oi -0.25 

43+ 41 



5+ ' 1 
U -°-0.! 



n 9+I.8 
z - z -1.5 



24011"! 

OtTQ + 140 



340+ 11(1 
J4U -100 

200+ 130 
zuu -120 



73.5/24 
46.0/22 

55.1/23 
52.6/23 
51.1/22 
19.0/21 

18.4/22 
12.2/21 



a This is the maximum temperature in the plT models 

b The normalisation of the plasma model in units of 10 _14 /(47rD 2 ) X Emission Measure 



of the components are assumed to have solar abundances, 
while the third has zero abundance (a bremsstrahlung spec- 
trum). The latter component acts as a means to dilute the 
strong line emission predicted by the other two temperature 
components. 

There is strongly significant additional low-ionization 
iron line emission. Including a narrow Gaussian emission line 
at 6.4 keV in the fit gives a reduction in xt to 1077/1157, 
with line equivalent width of ~ 140 eV. However, including 
the reflection continuum that is expected to accompany such 
line emission does not lead to a decrease in xl- The best 
fit reflected fraction is close to zero, with an upper limit of 
R < 0.7, too small to produce the observed 6.4 keV emission 
line. 

The low energy bandpass of ASCA means that these 
data are also sensitive to details of the absorption column 
along the line of sight, and physically some X-ray absorption 
should be present from the partially ionized wind seen in SS 
Cyg in outburst. The observed X-ray continuum spectrum 
can be described as a power law of photon index T ~ 2, so we 
use the the model absori in XSPEC, with the photon index 
fixed to this value. This leads to a highly significant decrease 
in xl to 1024/1154 for N H = 2.6x 10 21 cm" 2 at an ionization 
parameter of £ = 1.4 for a fixed temperature of 5 x 10 4 K. 
The underlying continuum form also changes, weighting the 
power law emissivity-temperature distribution even more 
towards the low temperature components, with a = 0.08, 
and the result is that the reflection component is now sig- 
nificantly detected, with R — 1.3. Models with the ionized 
absorption and 6.4 keV Gaussian line but without the re- 
flection continuum give x 2 = 1036/1155. Thus the reflection 
continuum is significantly detected in these data, but only 
when the intrinsic continuum form is properly modeled. 

The absorption in this partially ionized material has 
some effect on the derived bolometric luminosity of the X- 
ray emitting gas. Assuming that material cools down to the 
photospheric temperature of ~ 25 eV (Ponman et al 1995) 
gives a total flux from the hot plasma of 5.2 X 10 -11 ergs 
cm~ 2 s _1 (25 eV- 100 keV), while correcting for the ion- 
ized absorption increases this to 6.1 x 10 -11 cm~ 2 s _1 , and 
further correcting for the interstellar column gives a flux of 
11 x 10 -11 ergs cm -2 s -1 . This is almost identical to the 
(absorption corrected) hot plasma flux seen in the GINGA 



quiescent data integrated over the same band, showing that 
it is not necessarily true that there is less emission from the 
optically thin plasma in outburst, only that it is substan- 
tially softer. 

The spectral resolution and sensitivity of ASCA make 
it possible to investigate the abundances of the elements in 
the hot plasma emission separately rather than assuming 
that they all scale together. We replace the plasma emission 
model with one in which the abundance of each element 
is a free parameter (the cevmkl model in XPSEC, modi- 
fied as described in Section 4). The reflection abundance 
ratios are assumed to be solar, but the overall elemental 
abundance is a free parameter. This is fit separately from 
the plasma abundances as the plasma line emission may be 
distorted from that expected from coronal conditions (see 
section 5.1). Allowing the plasma iron abundance to be free 
relative to the other elements gives a significantly better fit 
{xl = 999.9/1152, for A Fc ~ 0.3, A rcst ~ 0.5). Progres- 
sively fitting for the abundance of Si and S (which are the 
next strongest expected lines) gives xt ~ 978.7/1151 and 
973.2/1150, for A Fc = 0.28, Ass = 0.72, As = 0.68 and 
A Icst = 0.41 (see table 4). A model in which all the other 
observable elements (O, Ne, Na, Mg, Al, Ar, Ca, Ni) are 
free gives xt = 964.2/1143 (C and N are fixed at solar as 
their Kq lines are below 0.6 keV so their abundance cannot 
be usefully constrained). Thus, from fitting coronal equilib- 
rium models (together with the reflection continuum and 
line and ionized absorption), the data require that the ele- 
mental abundances are not solar and that the abundances 
ratios are not solar either. The model with fewest parame- 
ters that can describe the data (where adding in extra el- 
emental abundance parameters does not significantly im- 
prove the fit), is one in which all the elements are tied to- 
gether except for Fe which is a free parameter, and Si and S, 
which are tied to each other. This gives xt = 973.4/1151, for 
Ap c ~ 0.3, Asi.s ~ 0.7 and A TCSt ~ 0.4. The spectrum was 
then truncated above 1.7 keV, so that the abundance of iron 
can be measured separately from the iron L complex. The 
contribution from reflection (continuum and line) can then 
be ignored, and the continuum parameters fcT max = 14.0, 
a = 0.15 and Asi,s = 0.70 were frozen at their best fitting 
value above. This gives xl = 244.8/221 for A Fc = 0.30 and 
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Table 3. Anomalous Outburst (ASCA) data 



model 


klQ (keV) 


A 


n£ 


a 


R 


EW (eV) 






xl 



IT 


o 4 o + o.io 


21 + - 05 


3.0 x 


10~ 2 












3799/1159 


2T 


flV+0 03 

u -°'-o.oi 
5.45 ±0.07 


0.51 ± 0.06 


3.4 x 
2.2 x 


10" 3 
10~ 2 












1331/1157 


3T 


"Sfi+ 04 
u -°°-0.09 

1.10 
6.64 


0.44 ± 0.06 


2.8 x 
3.0 x 
2.0 x 


10~ 3 
10~ 3 
10- 2 












1145/1155 


plT 


18.2 ±1.7 


n 4 o+0.06 
U - 48 -0.05 


3.1 x 


10~ 2 


0.47 ± 0.05 










1117/1158 


plT+g 


17.51™ 


n 44O.O6 
U ' 44 -0.05 


3.1 x 


10- 2 


0.45 ± 0.06 




140 ± 40 






1077/1157 


plT+R+g 




q 44+0-06 


3.2 x 


10~ 2 


46+ 06 

u -^ D -0.05 


O-OSlo.Sl 


140 ± 40 


2 6 +0 ' 7 
z - D -0.6 




1077/1156 


pia(plT+R+g) 


ia.e±i:S 


0.32 ± 0.05 


2.5 x 


10~ 2 


0.08 ±0.11 


1.3 ±0.7 


130 ± 40 


1.4 ±0.3 


1024/1154 


pia(plT+g) 


19 4+ 34 


0.41 ± 0.06 


2.4 x 


10- 2 


16+ 009 




160 ± 40 


2.1 ±0.6 


1 0+0.5 

1 -°-0.3 


1036/1155 



a This is the maximum temperature in the plT models 

b The normalisation of the plasma model in units of 10 _14 /(47rD 2 ) x Emission Measure 
c xlO 21 cm" 2 

d Defined as L/nr 2 where r is distance from material of density n illuminated by a source of luminosity L 



A Icst — 0.42, showing that the iron L complex is consistent 
with the iron K line strengths. 

The abundances in the reflector are harder to quan- 
tify, since they correlate with the amount of reflection (e.g. 
George & Fabian 1991). However, the observed 6.4 keV line 
emission is also dependent on the abundances in the reflec- 
tor and so provides an additional constraint. This is not 
straightforward, as the iron line equivalent width depends 
on the abundances of all the elements, not just iron, as the 
other elements provide opacity to the escaping iron line pho- 
tons (George & Fabian 1991). The functional dependence is 
also made more complicated by the fact that the refected 
line equivalent width also depends on the spectral shape of 
the illuminating continuum, as for hard spectra there is rel- 
atively more luminosity above the iron K edge that can be 
absorbed to produce the Ka line. However, knowing a line 
equivalent width from one particular illuminating spectrum 
we can scale this by the ratio of the absorped luminosity 
above the iron K edge to estimate the line equivalent width 
produced from reflection of any other spectral shape (Basko 
1978). The abundance effects can then be estimated from 
George & Fabian 1991, who show that for Morrison & Mc- 
Cammon (1983) abundance ratios then the line equivalent 
width scales roughly as for A a n < 1, and that a rela- 
tive over- or under-abundance of iron then gives a roughly 
proportional change in the line equivalent width. 

The iron line has an angle averaged equivalent width of 
110 eV with respect to a power law illuminating continuum 
of photon index T — 1.9 for Morrison & McCammon (1983) 
solar abundances (George & Fabian 1991). However, we are 
using Anders & Grevasse (1989) abundances, in which the 
iron abundance is 1.5 x larger than in Morrison & McCam- 
mon (1983), while all the other abundances are about the 
same. Thus with the Anders & Grevasse (1989) abundances 
we predict a line which is 1.5 x larger than this (Reynolds, 
Fabian & Inoue 1996), i.e. an equivalent width of 1.7 keV 
with respect to the reflected continuum. The ASCA spec- 
trum of SS Cyg has best fit spectral parameters of kT = 14 
keV, a = 0.15, which has 0.81 x fewer photons that can be 
absorbed by the iron K edge than the T = 1.9 power law, so 



predicting a line of 1.4 keV equivalent width with respect to 
the reflected continuum. If the abundance ratios in the re- 
flector scale together then for A a u ~ 0.3 the line equivalent 
width is reduced by a factor ~ 0.79, i.e. 1.1 keV. This is in 
strong contrast to the best fit model above (where the line 
and continuum from reflection are fit separately) in which 
the iron line equivalent width with respect to the reflec- 
tion continuum is ~ 0.6 keV. Fitting a reflection continuum 
in which the reflector abundances scale together with self 
consistent line emission leads to a strong increase in xl to 
979/1152 (i.e. A^ 2 of 6, compared to the fits in which the 
line normalisation is free) for an incident continuum in which 
the abundances in the hot plasma are all tied together ex- 
cept for Fe, and Si and S which are tied to each other. This is 
a strong indication that the iron abundance in the reflector 
is too high with respect to the abundance of the other ele- 
ments. The line is only consistent with that observed if the 
iron abundance in the reflector is subsolar in its ratio with 
respect to the other elements, and by an amount similar to 
that seen in the hot plasma line emission i.e. Ap c ~ 0.3 while 
A res t ~ 0.5. 

We replace the Anders & Grevasse (1989) 'solar' abun- 
dances with those of Morrison & McCammon (1983), which 
has a factor 1.5 lower iron abundance (see section 5.2), in 
both the plasma and reflection codes. Details of some of 
these fits are given in Table 5. We first allow the line to 
scale independently of the reflection continuum, and the re- 
flection spectrum abundances to scale independently of the 
hot plasma. With the abundance ratios tied at the (new) 
solar values we get xl = 1006/1153 with A ~ 0.4. There 
is a marginally significant change in the fit if the Fe abun- 
dance is fit separately, with xl = 1002/1152. However, this 
is driven by the relative overabundance of Si, and when Si 
is also free fit is significantly better with xl = 978.7/1151 
and the resulting Fe abundance is consistent with that of the 
rest of the elements. There is a further (marginally) signifi- 
cant decrease in xl to 973.6/1150 with S free to vary, but no 
other significant deviations from (new) solar element ratios. 
In particular, allowing the abundance of O to be free gives 
xl = 973.4/1149 for A Q = 0.42l°; 2 °. The model with fewest 
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Table 4. Abundances with respect to Anders & Grevasse (1989) 



kT 



(keV) 



A (plasma) 



A (reflector) 



R 



EW (eV) 



13.9. 



+6.2 



29+ - 07 
+0.15 
0.13 

S=0.65t°;2i 



Fe 

Si: 



=0.7CT 



rcst=0.40 ± 0.11 



0.15 ±0.11 



o.5: 



1.9 ±0.9 110 ±30 1.2 



+0.5 
0.3 



0.18 



+0.30 



973.2/1150 



a Maximum temperature in the plT models 
b xlO 21 cm" 2 

c Defined as L/nr 2 where r is distance from material of density n illuminated by a source of luminosity L 



parameters that describes the data is then one in which Si 
and S are tied together, but are free to vary from the Morri- 
son & McCammon (1983) solar ratios (xl = 974.9/1152, for 
Asi,s = 0.65, A rest = 0.37). The predicted iron line equiv- 
alent width with respect to the reflection spectrum is then 
~ 0.7 keV for kT = 14 and a = 0.15. This is a much better 
description of the fluorescent line seen in the data than that 
expected from the Anders & Grevasse (1989) abundances. 
We set the reflector abundances to be equal to A reat , and in- 
clude the self-consistent line emission in the reflection spec- 
trum, which gives a good fit at xt = 976.7/1154. 



4.4 All Data 

We fit a power law temperature-emissivity plasma model 
with Morrison & McCammon (1983) abundances, together 
with its reflection spectrum (also assuming these abun- 
dances) and a separate Gaussian line to all the data simul- 
taneously, with the abundances tied across all the datasets. 
An ionized absorber is included for the ASCA data, and 
the plasma emission and Gaussian line normalisations, the 
maximum temperature and temperature weighting index are 
free to vary between the spectra. This gives an excellent 
xt = 998.5/1197, showing that the elemental abundances 
are consistent across all the datasets, with Asi,s = 0.64 
and A rest = 0.40. The two outburst spectra are consistent 
with having the same amount of reflection, with no increase 
in xl when the two parameters are tied together, so that 
xl = 998.5/1198 for R = 2.2 in outburst and R = 0.7 in 
quiescence. This difference between the amount of reflection 
in outburst and quiescent is highly significant: tieing all the 
reflected fractions together leads to an increase in x 2 of 6.6. 

The best fit model spectra for the quiescent, GINGA 
decline from outburst and ASCA outburst data fit simul- 
taneously are shown in Figures 4, 5 and 6, respectively. 
The plasma parameters are a = 3.9, 0.8, 0.19, and kT max = 
13, 8.5, 12, respectively, again showing that most of the gas 
is hotter in quiescence than in outburst. 



5 DISCUSSION 

All the data strongly require that the line emission is sub- 
stantially weaker than expected from a coronal plasma with 
solar abundances. While this could simply indicate that the 
elements are subsolar, we first examine alternative possibil- 
ities. 



1> 



G 

o 

o 

£1 



0.5 1 2 5 10 20 

channel energy (keV) 

Figure 4. The best fit GINGA quiescent spectrum. 

5.1 Validity of the Coronal Approximation 

The data have been fit with coronal plasma models, which 
assume that radiative (and dielectronic) transitions are the 
dominant recombination processes balancing collisional ion- 
ization. All these processes involve only one electron in a 
free-bound (or bound-free) transition, so the resultant ion 
populations are independent of the electron density to first 
order. Most of the ions are in the ground state as the re- 
combination rates are generally very much larger than the 
collisional ionization rates. However, at high densities colli- 
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Table 5. Abundances with respect to Morrison & McCammon (1983) 



kT 



(keV) A (plasma) 



A (reflector) 



EW (eV) 



N, 



16.5 
13.9 



+5.6 
4.3 
+6.2 



0.42 
Fe=0.39 

Si=0. 

S=0.58 
rest=0.38 



+0.09 
0.06 
+0.09 
0.08 
,+0.18 
-0.16 
+0.18 
0.17 
+0.11 
0.10 



0.11 ±0.11 
0.17 ±0.11 



0.9 
0.5 



+4.1 
-0.5 
1.1 
0.25 



i.5: 



1.9 ±0.9 



0.7 



120 ± 40 
120 ± 40 



1.0 



± 0.6 

+0.6 



0.7 

0.36 



+0.4 
0.5 
0.0(1 
0.30 



1006/1153 
973.6/1150 



Maximum temperature in the plT models 
xlO 21 cm -2 

Defined as L/nr 2 where r is distance from material of density n illuminated by a source of luminosity L 
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Figure 5. The best fit GINGA decline from outburst spectrum. 



Figure 6. The best fit ASCA anomalous outburst spectrum 



sional excitation becomes important enough to rival radia- 
tive recombination, so that the ions are generally in excited 
rather than ground states. Collisional ionization out of the 
excited states, a process neglected in the assumption of a 
coronal plasma, can then modify the ionization balance. It 
also leads to suppression of the dielectronic recombination 
rate but this is not an important process at the tempera- 
tures considered here (e.g. Mewe 1990). For Hydrogen-like 
ions with atomic number Z in a plasma of temperature T 



(K), deviations from coronal ionization balance are less than 
a factor 2 when the density n < 4 x 10 4 Z 2 T 2 cm" 3 (e.g. 
Mewe 1990). Typical temperatures where the hydrogen-like 
state dominates for oxygen and iron are 10 7 and 10 s K, re- 
spectively, so these elements are out of coronal equilibrium 
for densities > 2 x 10 20 and > 2 x 10 23 cm" 3 . The lower 
limit on the density inferred for the boundary layer in the 
eclipse observations of HT Cas is > 2.5 x 10 13 cm" 3 (Mukai 
et al 1996). Thus there is no direct evidence against coronal 



© 0000 RAS, MNRAS 000, 000-000 



The X-ray Spectrum of SS Cyg 13 



equilibrium approximations for the ionization balance in the 
boundary layer in HT Cas in quiescence. However, the mass 
accretion rate in SS Cyg in quiescence is nearly two orders 
of magnitude larger than that in HT Cas, so its boundary 
layer density may well be substantially higher. The models 
of Popham and Narayan (1993) suggest a density of ~ 10 15 
cm~ 3 for conditions appropriate to the quiescent GINGA 
spectrum. Again, in outburst the boundary layer density is 
presumably higher still, but the complexity of the flow makes 
it difficult to estimate the density of the hard X-ray emit- 
ting material. Thus it seems unlikely that the material is as 
dense as required for coronal ionization balance to become 
invalid. 

A generally more restrictive assumption in the coro- 
nal models is that the material is optically thin. In HT 
Cas the eclipse measure can be used to set constraints 
on the quiescent emission properties. Assuming that the 
emission forms a corona of with outer radius (1 + S)Rwd 
then the column density through the boundary layer is 
N H ~ W 13 S 1/2 R WD < 3 x I0 21 cm~ 2 for the maximum 
value of S — 0.15. This is optically thin to electron scatter- 
ing and is also likely to be thin to all the H and He-like 
Lya transitions (Matt 1994; Mewe 1990). Even if these lines 
were optically thick, there is no viable line destruction mech- 
anism for these transitions (see e.g. the discussion in Done 
et al 1995) so the multiple scatterings have little effect. The 
plasma is possibly optically thick at iron L (Band et al 1990), 
but the fact that the iron abundance as measured from the 
L lines agrees with that derived from the K lines (see section 
4.3) argues against this strongly affecting the line emission. 

Self photo-ionization of the plasma is potentially impor- 
tant, since there is an intense X-ray flux. The effective ion- 
ization parameter for a diffuse source is £ = 4ncU /n where 
U ~ eSRwD/c, and e is the bremsstrahlung emissivity (see 
e.g. Done et al 1995). Thus £ ~ 1.8 x 10~ 14 tj ~ 18 for a 
density of ~ 10 15 and S — 0.15. This is not high enough 
to disturb the ionization balance set by the collisional pro- 
cesses, as the ionization states, even for Oxygen for a pure 
photoionized plasma are much less than those produced in 
a collisionally ionized plasma at kT ~ 2 — 10 keV (Kallman 
& McCray 1982). 

Another process that is neglected in the coronal models 
is three body recombination (the inverse process of colli- 
sional ionization). Here a free electron recombines with the 
ion, but the recombination energy is given to another free 
electron rather than being radiated as a free-bound tran- 
sition. Since this involves two thermal electrons, the rate 
depends on n 2 rather than n e as for the processes discussed 
above, so it becomes important only at high densities. Again 
it efficiently populates the upper excited levels of the ions, 
so that collisional ionization from excited states can be- 
come important and collisional excitation/de-excitation can 
change the expected line strengths. Cota (1987), quoted by 
Rees, Netzer and Ferland (1989), states that 3 body pro- 
cesses should be considered for all densities > 10 10 " 11 cm~ 3 , 
which is clearly within the range of densities which are im- 
portant here. However, it is unclear whether this can have 
an appreciable effect on the 2-1 line transitions, and Keady 
et al (1990) show that the coronal approximation for iron 
ion fractions in a 2 keV plasma holds up to densities ~ 10 20 
cm~ 3 . 

Coronal models also assume that the plasma is in equi- 



librium, and so do not apply where the physical conditions 
vary on timescales which are fast compared to the ion- 
ization/recombination timescales e.g. in young supernovae 
remnants (see e.g. Shull 1982). The collisional ionization 
timescale, ti „ = (nC) _1 , where C is the collisional ion- 
ization rate coefficient, is slowest for H-like iron (Shull and 
Van Steenberg 1982), where t io „ ~ 10 16 /(nT ' 5 ) at ~ 10 8 K, 
i.e. 0.1 sec at n = 10 13 cm~ 3 . This is much faster than the 
bremsstrahlung cooling timescale T coo i — 1.5 x 10 11 T ' 5 n _1 
sec for plasma at 10 K, and is even rapid compared to 
the observed X-ray variability (Yoshida et al 1992). This it 
seems likely that the plasma is hot and dense enough to be 
in ionization equilibrium. 

One density dependent process that will clearly affect 
the 2-1 line flux is collisional mixing of the 2s and 2p states. 
For an initially bare nucleus, a recombination cascade can 
end with an electron transition to the Is (1 2 S) level from 
either the 2s or 2p state (2 2 S or 2 2 P). The 2 2 S level cannot 
decay down to ground via a permitted electric dipole tran- 
sition as it violates the selection rule of AL — ±1. Its only 
decay path is a 2 photon electric dipole decay, leading to a 
continuum rather than an observable line. For recombining 
initially H-like ions the situation is somewhat more com- 
plex, as the cascading electron can have spin either aligned 
or misaligned with that of the electron in the Is state. Thus 
an excited 2s-ls He-like ion can be either 3 S, S P, 1 S or 1 P, 
of which the 1 S state decays via 2 photons. Thus the Ka 
line from H and He-like ions can be increased by collisional 
mixing out of the 2 2 S and 2 1 6' level, respectively, and this 
can strengthen the line by up to 50%, becoming effective for 
H-like ions of atomic number Z at densities greater than 
n ~ 1.8 x 10 4 Z 9 ' 3 (e.g. Netzer 1996). Thus it should be 
important for O at densities > 5 x 10 12 cm" 3 , and all ele- 
ments up to Si and perhaps S for densities of ~ 10 15 cm -3 . 
However, it is unlikely to affect iron as this requires densi- 
ties > 2 x 10 17 cm~ 3 . Thus in comparison to the plasma 
codes calculated with n = 1 cm we expect the observed 
Lya lines from iron in the hot plasma to be representative 
of the true abundance of iron, while those from the lower 
Z elements can be enhanced by this process, leading to an 
overestimate of their abundance by a factor of ~ 1.2 — 1.5. 

In short, while true coronal equilibrium conditions are 
not fulfilled in the boundary layer, none of the density, 
photo-ionization or optical depth effects discussed above are 
likely to cause the observed sub-solar abundance line inten- 
sities from the plasma. The one process that should be im- 
portant is enhancement of the Lya lines from H and He-like 
ions of low Z elements by a factor ~ 1.2 — 1.5 by collisional 
mixing transforming states that decay via 2 photons into 
Ka line. 

5.2 Abundances 

Since coronal models seem to give a good indication of the 
X-ray line equivalent widths, we can use these to estimate 
the abundances in SS Cyg. If collisional mixing affects ev- 
erything other than iron, then all the elements are a factor 
~ 3 sub-solar except Si and S, which are probably more like 
a factor ~ 2 sub-solar. However, this is not consistent with 
the intensity of the observed 6.4 keV fluorescent line from 
the reflector in the ASCA data, which is weaker in com- 
parison to the reflection continuum than is expected from 
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material in which all the elements are subsolar by a factor 
3. This strongly suggests that the relative abundance of iron 
with respect to the lighter elements is smaller than the 'so- 
lar' definition used here. A recent abundance compilation 
by Grevasse et al., (1996) gives Fe/H of ~ 3.2 x 1CP 5 com- 
pared to ~ 4.6 x 1CT 5 in Anders & Grevasse (1989), while 
most of the lighter elements have abundances within 10% 
of the Anders & Grevasse (1989) values. Ironically, these 
newest abundances are very similar to those of Morrison & 
McCammon (1983) for all the given elements. Without col- 
lisional mixing, the plasma in SS Cyg is then uniformly a 
factor ~ 0.4 below the 'solar' values of Morrison & McCam- 
mon (1983) except for Si and S, which are ~ 0.6x 'solar'. 



5.3 Amount of Reflection 

The determination of the amount of reflection is dependent 
on the assumed abundances. From the above discussion it 
seems that the abundance ratios in the plasma and reflector 
are closer to those of Morrison & McCammon (1983). Using 
these abundances for the reflector as well as the plasma gives 
the result for SS Cyg that R = 0.7 in quiescence, R — 2.2 
in outburst (see Section 4.4). The amount of reflection in 
quiescence is not consistent with the presence of an opti- 
cally thick inner disk (see Section 2.2 and 2.3). Either the 
inner disk is disrupted, or is optically thin. By contrast, the 
larger amount of reflection in outburst requires that there 
is additional reflecting material present compared to that in 
quiescence. This is consistent with a geometry in which the 
boundary layer extends over the whole white dwarf surface, 
illuminating both it and an inner disk. 



5.4 Ionized Wind 

A partially ionized absorber is seen in the ASCA outburst 
spectrum. The derived ionization state for £ ~ 0.2 — 0.5 
(tables 4 and 5) from the absori code (see section 4.0) is 
such that the CIV ionization fraction is at its maximum, 
with fciv ~ 0.1. This is the upper end of the range found for 
the blackbody boundary layer or disk illumination models 
of Hoare and Drew (1993), where fciv is typically 0.1 — 
0.001. The X-ray column density in the warm material in 
SS Cyg of N H ~ 10 21 cm" 2 gives M~2x 10" 12 M Q yr" 1 
using the spherical wind velocity law of Hoare and Drew 
(1993). Thus Mfciv ~2x 10 _13 M© yr" 1 , much smaller 
than that inferred from the UV resonance line modelling 
(e.g. Hoare and Drew 1993), perhaps because the ASCA 
data are taken towards the end of an outburst (Nousek et al 
1993), where the wind may be declining rapidly. We stress 
that the hard X-ray illumination does have important effects 
on the ionization structure of the wind, and that the X- 
ray absorption can give new constraints on the outflowing 
material which has hitherto only been extensively studied in 
the UV. 



5.5 Continuum Models 

The lack of robust theoretical models for the hard X-ray 
spectrum is a severe problem in interpreting the data. The 
calculations of Narayan and Popham (1993) are the only 



numerical results for the quiescent spectra that can be rig- 
orously tested. However, these are restricted by the 1 di- 
mensional approximations which imply vertical averaging of 
what is most likely to be a plasma with a stratified tem- 
perature structure. In these models the expected spectrum 
is very like that produced from material that is impulsively 
heated which then cools via bremsstrahlung alone in a con- 
stant pressure/gravitational field. This spectrum is incom- 
patible with the GINGA quiescent spectrum, in which the 
expected cooling gas is not present. Either these models for 
the quiescent spectrum are wrong, and there is some addi- 
tional process heating of the gas, or the cooler components 
are masked by absorption. The ionized wind is is unlikely to 
be important in quiescence, so the only potential absorber is 
an optically thin inner disk, which could attenuate the emis- 
sion from the observable part of the hard X-ray boundary 
layer below the equatorial plane. However, this absorption is 
not required by the data - an equally acceptable description 
of the spectrum is an X-ray source where the temperature- 
emissivity relation is weighted towards higher temperatures 
than expected from simple cooling models. 

In outburst there are no detailed models for the 
temperature-emissivity structure that can currently be com- 
pared with the data. However, the integrated luminosity 
from the hot gas is similar in both outburst and quiescence, 
although its fractional contribution to the bolometric lumi- 
nosity is much lower in outburst. The maximum temper- 
ature component of the emission is similar in both quies- 
cence and outburst at ~ 10 — 15 keV (see section 4.4). As 
expected, these are always less than the maximum tempera- 
ture of ~ 20 keV produced by a strong shock from material 
accreting onto a white dwarf mass of O.8M and radius of 
7 x 10 8 cm (Frank et al 1992). 



6 CONCLUSIONS 

We find a significant reflection signature in all the SS Cyg 
spectra. This is stronger relative to the illuminating plasma 
emission in outburst than in quiescence, consistent with 
models in which the inner disk is optically thick in outburst 
but optically thin or non-existent in quiescence. The amount 
of reflection in both quiescence and outburst also supports 
the idea that the hard X-ray plasma forms a corona over 
the white dwarf surface, rather than being confined to an 
equatorial band where the disk/star interaction takes place. 

The ASCA outburst spectrum shows significant absorp- 
tion from partially ionized material, which is probably the 
outflowing wind inferred from the UV resonance line profiles. 
The X-ray absorption gives an additional diagnostic win- 
dow on the wind properties, and suggests rather low mass 
loss rates of ~ 10~ 12 M Q yr _1 . The hard X-ray emission 
strongly affects the ionization structure of the wind, giv- 
ing much larger CIV ion fractions than previous blackbody 
boundary layer photoionisation models. 

The ASCA data also give a good determination of the 
amount of line emission from the plasma, which are a factor 
~ 2.5 weaker than expected for solar abundances. We crit- 
ically examine the assumptions of coronal equilibrium but 
conclude that the line strengths are representative of the 
abundances. 

The intrinsic plasma emission, when separated from the 
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effects of reflection and (in outburst) ionized absorption, can 
then be described by hot plasma with a continuous (power 
law) temperature-emissivity distribution. In quiescence this 
distribution is strongly weighted towards the highest tem- 
peratures, with very little cooler material present, so that 
the emission can in practice be modeled by a single tem- 
perature plasma. The lack of cool components is in conflict 
with the predictions of the advective boundary layer models 
of Narayan & Popham (1993), and with simple models of 
cooling gas! One possible way to reinstate these models is 
for the quiescent spectrum to be distorted by partial absorp- 
tion, such as would be produced if the inner disk is optically 
thin to electron scattering but optically thick to photoelec- 
tric absorption. The flux from the boundary layer below the 
disk is then absorbed, while that from above the disk is not. 

The outburst spectrum by contrast is much softer. It is 
dominated by the cool components and cannot be modeled 
as single temperature emission, but the total luminosity of 
the X-ray emitting gas is very similar to that in quiescence. 
We review current theoretical models for the hard X-rays in 
outburst, but none of these address the probable difference 
in the physical nature of the viscosity between the disk and 
the boundary layer, or the inherently 3 dimensional nature of 
the problem. We urge further theoretical modelling, giving 
the predicted, spatially integrated temperature-emissivity 
distribution, so that a better understanding of the spectra 
of dwarf novae can be developed. 
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8 APPENDIX 

The angle r) between the X-ray emitting source and the line 
of sight is then given by 

cos rj = sin i sin 9 sin + cos i cos 9 



an inner disk or the extent of the boundary layer. Thus the 
reflection component from the white dwarf surface in SS 
Cyg is approximately that expected from an isotropically 
illuminated slab viewed at 60°. 



The reflection spectrum from each small segment of the 
white dwarf underlying the boundary layer can be approx- 
imated by an isotropically illuminated slab viewed at incli- 
nation rj. Thus the average inclination between the reflector 
and the observer is given by averaging ij over all angles at 
which the boundary layer emission can be seen. If the disk is 
opticaly thick and extends down to the white dwarf surface 
(figure 2a) then this is given by 



2 



tt/2 fir/2 



v(i,P) = T7 / # / d9 V (9,cb,i) sin 9 

1V Jo Jtt/2-P 

+A f d4> ( d0tj(M>») sm6 (1) 

J-4>min Jtv/2-0 

where the normalisation N is given by equation (1) but with 

all the r)(9, <f>, i) terms removed from the integrand, 

and 



= Tv/2...(3>i 
"tan P 



tan i 



■ P<i 



arctan 



( - 1 ) 

\sin0tani/ 



This needs to be integrated numerically. The inclination 
angle of SS Cyg is thought to be i — 37° (Cowley et al 
1980). For the two limiting values of /3 -> 0° and /3 -> 90° 
the angle of the emission to the line of sight is 67° and 53° , 
respectively. 

Alternatively, if the inner disk does not obscure the 
boundary layer (figure 2b) then the mean inclination angle 
is 

IT/2 /■7T/2 + /3 



«L d4 L 



r)(i,P) = - I d( t> I Mr >(8, 4>,i) sine 

/2-fl 



+ Tf d(t> I d9ri{9,4>,i) sin 9 

N I" Jtt/2-P 



+ 1 / d4» I d9 V (9,4>,i) sin 9 (2) 

JV ' - Itv/2-P 



where S ma x and m i n are given as before, and the normalisa- 
tion N is given by equation (2) with all the rj(9, <j), i) terms 
removed from the integrand. 

Again, the limiting case for (3 — > is a mean inclination 
of 53° (the presence or absence of a disk makes no difference 
if the boundary layer is limited to the equatorial plane), 
while that for [3 -> 90° is 57°. 

Given that all these limits are very close to each other, 
the mean angle at which we see the boundary layer in SS Cyg 
is approximately 60°, irrespective of whether or not there is 
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